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Abstract. New results are presented for the stress conjugate to arbitrary Eulerian strain mea- 
sures. The conjugate stress depends on two arbitrary quantities: the strain measure /(V) and 
the corotational rate defined by the spin ft. It is shown that for every choice of / there is 
a unique spin, called the f-spin, which makes the conjugate stress as close as possible to the 
Cauchy stress. The f-spin reduces to the logarithmic spin when the strain measure is the Hencky 
strain In V. The formulation and the results emphasize the similarities in form of the Eulerian 
and Lagrangian stresses conjugate to the strains /(V) and /(U), respectively. Many of the 
results involve the solution to the equation AX — XA — Y, which is presented in a succinct 
format. 



The notion of stress and strain are interhnked, regardless of the existence of a strain energy 
function. At the most basic level they are related by mechanical power, the rate of work per 
unit current volume of material. 



Here cr is the Cauchy stress and D the stretching tensor. This work-conjugate relation is 
independent of any notion of a reference configuration, although it is useful to introduce one. 
Let F be the deformation gradient between the current and reference states, and let T and E 
be the stress and strain associated with the reference state. T and E are mutually conjugate if 
they satisfy 



where the factor det F arises from the change in volume between the current and reference 
descriptions. In fact, eq. (jl-2p is usually taken as the starting point for determining stress. The 
choice of the strain E is not unique, but once chosen it fixes the definition of T through the 
work conjugacy (jl-2p . It is strange but true that the same simple connection does not apply 
to the relation between current or Eulerian strain and the Cauchy stress. The difficulty is in 
the definition of strain, say e. What e is such that e = D? It turns out that this question is 
incomplete and that we must broaden it and seek the strain for which e= D, where ° signifies a 
corotational rate. Actually, the corotational rate itself also has to be found. Fortunately, both 
the strain and the rate have been determined: Xiao et al. pO] showed that the unique solution 
is obtained by the Hencky strain InV in combination with the logarithmic rate. But we are 
getting ahead of ourselves. 

It is evident that work-conjugacy is simpler for reference or Lagrangian stress and strain than 
for their counterparts in the current or Eulerian configuration. Note that the distinction between 
Lagrangian and Eulerian is made explicit by the polar decomposition F = RU = VR: quantities 
associated with or defined by U and V will be called Lagrangian and Eulerian, respectively. 



1. Introduction 




(1-1) 




(1-2) 
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It is instructive to review work-conjugacy for Lagrangian stress and strain. The starting point 
is the fact that the stretching tensor D is the symmetric part of FF~^. Let the strain be chosen, 
quite generahy, as E = /(U) where the function / is sufficiently smooth, then eqs. (Il-ip and 
(fr2l) imply 

tr (T [V/(U)] U) = tr(crD) det F. (1-3) 
The gradient V/(U) is a fourth order tensor function which will be described later. At the same 
time the kinematic quantities, strain rate U and stretching D, may be related quite easily (see 
Appendix lA]) 

U = 2(UKI + IKU)-^(UK1U)R*DR. (1-4) 

Using the independence of D, eqs. (jl-3p and (|l-4p imply, formally at least, that the stress 
conjugate to the Lagrangian strain /(U) is 

T=(V/(U))-'t«, (1-5) 

where T^^\ sometimes called the Biot stress or the Jaumann stress, and S, the second Piola- 
Kirchhoff stress tensor, are 

T(^) = ^(UKI + IKU)S, S = F"VF"*detF. (1-6) 

We have used the symmetry of T and certain commutative properties to express the stress in 
()l-5p as a fourth order tensor acting on T*-^^. The tensor product notation, Kl, explained in the 
next Section, is used throughout as we find it makes results more transparent. Equations (jl-5p 
and (jl-6p embody work-conjugacy for arbitrary Lagrangian strain /(U). 

Although the notation in eq. (jl-5p might be unfamiliar the result is not, see jl3( eq. (3.5.31)]. 
The fourth order gradient tensor V/(U) is discussed in detail by Norris [12]. In particular, it 
is positive definite, symmetric and invertible for any strain measure function [6]. Examples will 
be presented for the Seth-Hill strain measure functions. 

For instance, the stress T^™^ associated with /(™)(U) is 

t(°) = ( V In U) t(i) = d X K U^-^ t(^) , (l-7a) 

Jo 

T(^) = -(U^/3 K I + U^/^ M U^/^ + im U^/^) t(^), (l-7b) 
3 

T(5) = i(VUKI + IK\/U)T(^), (l-7c) 

t(2) = S, (l-7d) 
r^i-m) ^ (jjm ^ u*") T^*") . (l-7e) 

Some of the conjugate stresses listed are well known, e.g. m = 1, 2, —2 [6l[13], and (jl-7ep follows 
from p. 158]. Identities p-7ap - ()l-7cp and the others will become evident later. The second 
identity in ()l-7ap follows from [12]. We note that the Piola-Kirchhoff stress is conjugate to 
E = ^(U^ — 1), the Green strain, which is typically used in applications. 

A principal objective of this paper is to find analogous expressions for the Eulerian stress r 
conjugate to the strain e = /(V) where V is the right stretch tensor and the function / is again 
arbitrary. We also require that the Cauchy stress be included among the Eulerian stresses, just 
as the Piola-Kirchhoff stress appears naturally for the Green strain. 
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Unlike the Lagrangian strains, /(V) is not an objective tensor [2j, and it is known that this 
strain rate does not, in general, possess a conjugate stress [lOj. This difficulty can be avoided by 
defining conjugacy in terms of corotational strain rates. The corotational rate of a symmetric 
second order tensor A(t) is 

A= A + Afi - OA, 

where the skew symmetric tensor fl is called the spin. Xiao et al. [22^ [21] showed that an 
objective spin has the general form 

r2 = W+P(V)D, (1-8) 

where W is the skew symmetric part of FF^^ and P is an isotropic fourth order tensor-valued 
function of V. Lehmann and Liang [8j showed that using the rate associated with R, i.e a 
corotational spin equal to the "twirl" = RR*, the Eulerian and Lagrangian stresses conjugate 
to /(V) and /(U) are related by r = RTR*. This relationship simply rotates the Lagrangian 
stress, but does not reproduce the Cauchy stress for any choice of /. 

The fundamental relation for Eulerian conjugate stress is based on the finding of Xiao et al. 
[m that 

°log 

ln(V)= D, (1-9) 
where oj^g denotes an objective corotational rate defined by the logarithmic spin 

[19]. We 

will discuss in detail, providing a new derivation and representation, and comparison with 
ri^. The relationship (jl-9p allows us to define a class of work-conjugate Eulerian stress-strain 
pairs for all /(V) that includes the Cauchy stress. However, it should be borne in mind that 
the logarithmic rate is but one from a continuum of possibilities. 

A second objective of this paper is a generalization of eq. (jl-9p to arbitrary strain measure 
/(V). Xiao et al. [T9] proved that D is recovered only from the Hencky strain ln(V) combined 
with the logarithmic spin; no other strain measure can yield D no matter what spin is used. 
Here we will show that for a given strain measure /(V) there is a unique spin which provides 
the best approximation to D, and the corresponding conjugate stress is the best approximation 
to the Cauchy stress. 

lA. Summary of principal results. Our first main result is: 
Theorem 1. The stress conjugate to the Eulerian strain /(V) is 

T=(V/(V))-'r«(n), (1-10) 

where r^^^ depends on the corotational rate used, 

T^'^\ft) = [V^ K I -M S (V^ K I - I S V^) P(V)] K I I K V) cr. 

We will explain this result in detail, and provide alternative representations for eq. (ll-lOp . 
We note at this stage the similarity in form between the relations (11-51) and (ll-lOp . In particular, 



4 



ANDREW N. NORMS 



the stress r^™) conjugate to the strain /("-^CV) is 

r(°) = (VlnV)"VW = / dxV^ K V^-^rW, (1-lla) 

Jo 

r(i) = i(v2/3 M I + V^/^^ M V^/^ + j ^ v^/^) r^^), (1-llb) 



3 



r(2) = 2(VKI + IK V)~V(^), (1-lld) 

t(-™) = (V" K V") t("). (1-lle) 

The second principal result defines a new spin defined by the Eulerian strain measure. 

Theorem 2. For every Eulerian strain measure /(V) there is a unique corotational rate which 
minimizes the difference between the conjugate stress and the Cauchy stress. The rate is defined 
by the f-spin ftf = W + P-^D which depends upon the function f via the fourth order projection 
tensor 

= (VKI - IK V)*[(V/(V))"^ - (VKI + IK V)-1(V2 KI + IK V^)], 

and A* denotes the pseudo-inverse (or Moore-Penrose inverse) of the tensor A. The conjugate 
stress using the f-spin is 



where Aj are the principal stretches, Vj the principal dyads, i.e. the eigenvalues and eigentensors 
of\, and the eigen-index n G {1, 2, 3} is the number of distinct eigenvalues. The conjugate stress 
is minimal in the sense that \t — a\ > \cTf — a\ for any other corotational rate. 

The pseudo-inverse is a unique quantity and will be defined in detail later. 

The logarithmic spin |20j is a very special case of the f-spin. It is clear from eq. ()l-12p that 
when f{x) = Inx and the f-spin is used then the conjugate stress is simply the Cauchy stress, 
i.e. ct'"^ = <T. Note that cr is recovered as t^^\ the stress conjugate to the Hencky strain InV. 
No other spin reproduces the Cauchy stress as the conjugate of any strain [19]. This emphasizes 
the mutual relation between the Hencky strain and the logarithmic spin. 

IB. Review and plan of the paper. No attempt is made to summarize the considerable 
literature on work-conjugacy, strain measures and associated stresses, although two introductory 
reviews are worthy of mention. Curnier and Rakotomanana [T] provide an instructive overview 
of strain measures and conjugate stresses, with extensive references to the literature prior to 
1990. A more concise but in-depth description of work conjugacy and its implications is given 
by Ogden [13]. These reviews and most of the work prior to 1991 dealt with stress conjugate to 
Lagrangian strain measures, although there had been some relevant work on quantities related 
to InV. For instance, Fitzgerald ^ considered the stress conjugate to InV in the context of 
hyperelasticity. Hoger [7] derived expressions for the rate of change of U, which subsequently 
proved useful for Lehmann et al. [9l [8] when they considered V specifically. The focus here is 
on Eulerian strain and its work-conjugate stress, and builds upon developments in the 1990s. 
Lehmann and Liang [8j introduced a clear procedure to extend the idea of work-conjugacy 
to strains whose rates are not objective in a fixed frame, see also [9]. The idea, reviewed in 
Section m permits the use of corotational rates. This is especially important for Eulerian strain 
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measures since Xiao et al. |19] proved that the only way to obtain D with Eulerian strain is as 
the logarithmic rate of In V. We emphasize that the pointwise rate of working ib is the focus 
as we consider its implications for pointwise stress based on different definitions of strain. No 
assumptions of material homogeneity, isotropy, or otherwise is assumed or required. 

Xiao, Bruhns and Meyers provide the most complete and thorough analysis of Eulerian con- 
jugate stress and strain. In a series of groundbreaking papers [El [22l [2lj culminating in |2Q] 
these authors showed that the notion of conjugate stress is just as relevant to Eulerian strain 
as it is for Lagrangian strain. Because the role of the rate, or spin, is central to the Eulerian 
problem but is absent from Lagrangian work-conjugacy, it is essential to have a thorough un- 
derstanding of the possible spin tensors and their dependence on quantities such as D, W and 
V. Once this is understood then the form of the conjugate stress becomes apparent. Xiao et 
al. [20] derived expressions for the Eulerian conjugate stress for arbitrary strain measures / and 
for arbitrary permissible corotational strain rate. Their subsequent work has highlighted the 
role of the logarithmic rate and the Hencky strain in applications to hyper elasticity and other 
constitutive theories, see jl8j for a thorough review. 

This paper presents new results which extend the work of Xiao et al. in several directions. 
The introduction and discovery of the role of the f-spin shows that there is a certain unique 
conjugate stress associated with every Eulerian strain measure. The dual formulation for the 
Eulerian and Lagrangian conjugate stresses in eqs. ()l-5p and (jl-lOp further emphasizes the 
similarities in the two descriptions. The formulation throughout is in direct tensor notation, 
which we believe makes the results more transparent. 

The plan of the paper is as follows. The notation is introduced in Section [21 where the gradient 
and the pseudo-inverse of a tensor are defined. Corotational strain rates are discussed in Section 
[3] and some basic results for Eulerian strain measures are derived. The f-spin is introduced and 
discussed in Section [H It is shown that the corotational rate defined by the f-spin, or f-rate, has 
certain unique and desirable properties. The main results for conjugate stress-strain pairs are 
deduced in Section [5l 



2A. Preliminaries. We will be dealing with tensors of second and fourth order. Second order 
tensors act on vectors in a three dimensional inner product space, x Ax with transpose A* 
such that y • Ax = x • A*y. Spaces of symmetric and skew-symmetric tensors are distinguished, 
Lin = Sym © Skw where A G Sym (Skw) iff A* = A (A* = —A). The inner product on Lin is 
defined by A • B = tr(AB*). Products AB GLin are defined by y • ABx = (A*y) • Bx. 

Psym is the space of positive definite second order tensors. When dealing with functions of 
a symmetric tensor it is often useful to rephrase the functional form in terms of the spectral 
decomposition: 



2. Tensors functions and the pseudo-inverse. 




(2-1) 



where Aj GPsym and n < 3 is the eigen- index. Thus, 



n 



/(A) = j;/(a,)A,. 



i=l 

The Poisson bracket of two second order tensors is 



{A,B} = AB -BA. 
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Lin is the space of fourth order tensors acting on Lin, X — > AX with transpose A* such that 
Y • AX = X • A*Y for all X, Y GLin. The vector space may be decomposed Lin = Sym® Skw 
where Sym and Skw denote the spaces of symmetric ( A* = A) and skew-symmetric (A* = — A) 
tensors, respectively. Any A G Lin can be uniquely partitioned into symmetric and skew parts: 
A = A(+) + A(-), where A(=^) = ( A ± A*)/2. The identity I satisfies IX = X for all X G Lin. 
The product AB G Lin is defined by Y • ABX = ( A*Y) • BX. Psym is the space of positive 
definite fourth order tensors: A G Psym iff X • AX > 0, for all nonzero X G Sym. 

The square tensor product X Kl Y, LinxLin^ Lin, is defined by [14j 

(XKY)Z = XZY*, VZGLin. 
In particular, we note the property 

(A Kl B)(X K Y) = (AX) M (BY). 

2A.1. The tensor gradient function and its inverse. The gradient of a tensor function /(A) is a 
fourth order tensor V/ G Lin defined by 

V/(A) X = lim 1 [/(A + eX) - /(A)] . (2-2) 

We make extensive use of the following representation which uses the spectral form of A, 

v/(A)=f: :^m^/Ma.^a„ 



where the ratio becomes f'{ai) for i = j. The formula (j2-2p for the first derivative is well known, 
e.g. |13tll7|. Norris [12j provides formulas for the n*'* derivative of a tensor valued-function. We 
define the inverse tensor function A/(A) G Lin by 

n 

where the ratio is l//'(aj) for i = j. The definition of A/(A) is problematic if f'{ai) vanishes, 
but we preclude this possibility next by restricting consideration to strictly monotonic functions: 
strain measure functions. 

2B. Strain measure functions. The function / is a strain measure [6l [15] if it is a smooth 
function / : — > M which satisfies 

/(I) = 0, /'(I) = 1, /' > 0. 

It may be shown |12j that the gradient of a strain measure function and its inverse are positive 
definite fourth order tensors, i.e. V/(A), A/(A) G Psym. We restrict attention to strain 
measure functions for the remainder of the paper. 

2C. The pseudo-inverse. For A G Psym consider the equation 

{A,X} = Y, (2-3) 

for the unknown X in terms of Y. It is assumed that Y is either symmetric or skew and X is 
of the opposite parity [3J. The equation can be written AX — XA = Y, or 

JI(A)X = Y where J(A) = AK1I-IKA. 
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We will only consider JI(A) for symmetric A, implying J G Sym and J maps Sym^Skw and 
Skw^Sym. Therefore, J does not possess eigenvalues, eigenvectors or an inverse in the usual 
sense. 

The unique solution of the tensorial equation (j2-3p is [12] 



X= r(A)Y = (AKI-IK A)*Y. (2-4) 
The pseudo-inverse, or equivalently the Moore-Penrose inverse, J*, is defined such that 

jirji=ji, rir = r. (2-5) 

The spectral forms of JI(A) and its pseudo-inverse are 

n n 

JJ(A) = ^{ai- aj)A, M Aj, J* (A) = ^ (a^ - aj)-^Ai M Aj. 

i,j=l i,j=l 

These clearly satisfy ()2-5p . 

Further insight into the pseudo-inverse is gained by introducing the set of < 6 fourth order 
tensors associated with A GSym, 

(AjMAj, I = l,...,n, , , 

A/ = < ' ' ' > 2-6 

[AiMAj+AjMAi, I = n + l,...N. ^ ^ 

A'" = 6 for n = 3 and the indices / = 4, 5, 6 correspond to = (2, 3), (3, 1), (1, 2), respectively. 
Similarly = 3 if n = 2 and = 1 if n = 1. Note that 

I = |;A,, A,Aj = |q'^' ^j^j (2-7) 

The identity I = I Kl I implies the partition of unity (|2-7P i . and it may be readily checked that 
the A/ satisfy the orthogonality conditions (|2-7p 9. 
The pseudo-inverse satisfies 

n N 

JI*JI= jr = I-^ A/ = ^ A/. (2-8) 

/=1 I=n+1 

This is never equal to the identity I, which is the property that distinguishes the pseudo-inverse 
from the standard notion of inverse. Further properties of the pseudo-inverse are presented in 
Norris [12]. The explicit solution of eq. (12-41) can be expressed in a variety of ways without the 
use of fourth order tensors. Perhaps the simplest is the recently discovered solution of Dui et al. 

X= (3A'^ - i(trA'^)l)~\2A'Y + YA'), 
where A' is the deviatoric part of A. 



3. Kinematics 



3A. Basics. The polar decomposition of the deformation gradient is F = RU = VR where 
R gS0(3) satisfies RR* = R*R = I and the right and left stretch tensors U and V are positive 
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definite and related by V = (RKIR)U. Tlie fundamental Eulerian strain can be taken as either 
V or its square B = = FF*. The spectral representations of V and B are 

n n r . _ . 

V = J2\^^, B = J^ftV„ V,V, = ^*' (3-1) 

i=l 1=1 I ' ^ 7^ J) 

where Aj > and f3i = Xf. 

The rate of change of B is B = LB + BL* where L = FF"^ Let D eSym and W eSkw be 
the symmetric and skew-symmetric parts of L, respectively. Thus, L = D + W and B can be 
expressed 

B = (IK1B-BKI)W + (IKB + BK1I)D. (3-2) 
We will find this form useful for deriving more general strain rates. 

3B. Co-rotational rates. Let A(t) be a symmetric second order tensor, and ft is skew and 
arbitrary. Define the corotational rate 

A=A + {A,n}, fieSkw. (3-3) 
For any Q{t) ESkw we can identify a rotation Q(t) gS0(3) such that 



QAQ* = Q A Q*. 

Differentiating the left member and using (j3-3p for the right member implies that CI = — Q*Q. 
Hence, Q must satisfy Q = — Qfi, with solution unique up to a rigid body rotation. The coro- 
tational rate may therefore be interpreted as the Lie derivative with respect to spatial rotation 
defined by Q(t). Thus, let (p define the mapping (rotation) x Qx, then the corotational rate 

oj 

The Jaumann rate A defined by 12 = W corresponds to P = in eq. (jl-8p . Using the latter 
formula to parameterize the spin Q allows us to express the general corotational of A as 

A=A +(AKI-IK1 A)P(V)D. (3-4) 

Equation (l3-2p implies that the Jaumann rate of B is B= (I Kl B + B Kl I)D. The general rate 

o o o o 

B then follows from (j3-4p . and V can be determined from the identity B= (V Kl I + I Kl V) V. 
In summary, the general form of the corotational rate of the fundamental Eulerian strains are 

B=[BKI + IK1B + (BKI-IKB) P(V)] D, (3-5a) 
V = [(V Kl I + I K V)-^{V^ K I + I Kl V^) + (V K I - I K V) P(V)] D. (3-5b) 

3C. Spins. Many candidates have been considered from the infinity of possible spins [2]. For 
instance, the polar spin 

ft^ = RR*, (3-6) 

corresponding to Q = R*, is useful as a comparison spin. Other common spins [22] are Ct^ 
defined by the "twirl" of the Eulerian principal axes and fl^ related to the Lagrangian principal 
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axes. It is shown in Appendix O that 11" = W + P"D, a = R,E, L, where 

= (IK V- VKI)(ISV + VKI)"\ (3-7a) 

P^ = (I K - K I)*(I K + K I), (3- 7b) 

P^ = (I S - S I)* 2V M V. (3- 7c) 

The three spins fi^, fl^ and fl^ are related by fl^ - = fl^ - W, see Appendix S The 
fourth order projection tensors are therefore connected by P^ — P^ = p-^, and we note the 
additional relation P^ + P^ = P^ , which is readily verified. 

The most general form of the isotropic tensor- valued function P € Sym involves three isotropic 
scalar functions i^i , 1^2 , [22j , 

P(V) = (V K I - I K V) [z^i I + zy2(V K I + I K V) + i^sV m V] 

n 

= I^KjV,KVj, (3-8) 

where 

Pij = (Aj - Xj) [1^1 + (Aj + Xj)u2 + XiXji^-s] , t^k = T^k{Ii,h, h), k = 1,2, 3. 

Here, Ii, I2, 13 are the invariants of V: Ii = tr V, I2 = ^/f — ^ tr V^, I3 = det V. 
The corotational rate of V can now be written 

V= QD, (3-9) 
where the fourth order tensor Q G Sym follows from ()3-5bp . 

Q = (VKI + IK V)-i(V2 K I + I K V^) + (V K I - I K V) P(V) 

n 

and 

A^ + A^ 

Qij = (A. - Xj)p^j + (3-10) 

+ Aj 

3D. Eulerian strain measures. The Lagrangian Seth-Hill strain E^"^^ = m~^(U"^ — I) corre- 
sponds to f{x) = /^""^(x). We define the analogous Eulerian strain 

e(™) =/(™)(V) =m"i(V™-I), (3-11) 

and note in particular the Hencky strain e^''^ = InV. Other examples include 

e{i)=V-I, e(2) = l(B-I), e(-i)=I-V-\ e^-^) = 1(1 - B^^). 

3E. Eulerian strain rates. We now present some identities for the corotational rates of Euler- 
ian strains. These will prove useful later in deriving conjugate Eulerian stresses. The first 
identity applies to arbitrary strain measures: 

Lemma 1. The corotational rate of any Eulerian strain measure f(V) is 

/(V)=[V/(V)] Q(V)D. 
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The proof is a simple application of the chain rule using eq. (j3-9p for V. This separates the 
dependence on the strain measure / from the dependence on the particular corotational rate 
used, which determines Q. 

The second identity connects the strain rate with the Hencky strain: 

Lemma 2. The strain rate of any Eulerian strain measure /(V) can he expressed in terms of 
the Hencky strain rate as 

/(V)=[V/(V)](AlnV) InV. 

The proof is a straightforward generalization of the chain rule of differentiation |22l Thm. 2] . 
Let M = In V and /(V) = /(M) then, 

/(M)= Vm/(M) M= V/(V) (VmV) M . 
But the fourth order tensor VjyjV is just the inverse of since V-y"V = I. 

4. The f-spin and the logarithmic spin 

4A. Strain rate and the stretching tensor. In order to make the connection between the 
kinematics and the power w we must relate some strain rate to the stretching tensor D. A 
general connection can be found by starting with the rate of change of an arbitrary tensor 
valued function of B. Thus, 



n 



/(B) = [V/(B)]B = ^ 



«J=1 



where the temporary definition f{x) = /(x^) is used, so that /(V) = /(B). Reverting to /(V) 
and using Pi = X'^, the rate of change of the associated function of V is 

/(V) = E V, m V,B, (4-1) 

i,j=l * J 

where the ratio becomes /'(Aj)/(2Aj) for i = j. Substituting B into (j4-ip and using the filtering 
properties of Vj, such as (Vj M\'j){IMB) = PjYi M Yj, gives 

/(V) = E (/(A,) - f{x^))v, m v,w + E ^^^'^ I i^^'^ (A + /3,)V. M V, D 

• • 1 ■ ■ T Pi Hi 

1,3=1 «J=1 

n ^2 _|_ ^2 

= { w, /(V)} + E 3^h^(^(^^) - ^(^^■)) ^ ^ 

i,j=l * 

Adding and subtracting terms, this becomes 

" + 1 

/(V) = D + {W, /(V)} + E [-^^ - f^^^^ _ f^^^^ ] {fi>^^) - /(A.)) V. M V, D, (4-2) 
where D is a modified version of the stretching tensor, 

n 

D = D + E [Ai/'(A.) Vi D. (4-3) 
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Note that the double sum in ()4-2p excludes the i = j terms. We can therefore rewrite it in a 
form suggestive of a new corotational rate, 

/(V) = D + {n-^,/(V)}, ftf = W+FfB, (4-4) 

where GSkw is called the f-spin, and its fourth order projection tensor is 

= (V K I - I K V)* [A/(V) - (V S I + I K V)-\V^ MI + IM V^)] . (4-5) 

Alternatively, P-^ can be expressed in the form (|3-8p with matrix elements 

/_ 1 Af + A^ 

Note that fi-^ can blow up, but the action [O,^ ,g(V)^ remains finite for any differentiable 
function g, including /. In particular, the f-spin is an objective material spin in the sense 
defined by Xiao et al. [22j . 

°j 

The corotational rate associated with the f-spin is defined in the usual manner as g{\) = 
g{y) + {g(V), fi-^}. The reason for introducing this new rate is 

/(V)= D, 

which follows from eq. (j4-4p . This shows that for a particular choice of spin the corotational 
rate of an arbitrary strain measure /(V) is related to the modified stretching tensor D. The 
important point is that this is the closest, in a sense to be defined, the strain rate can get to the 
actual stretching tensor D. These ideas are made concrete through: 

Lemma 3. For any objective corotational rate 

I /(V) -Dp = I /(V) -Dp + |D - Dp, (4-6) 
where D is the modified stretching tensor defined by eq. ^-3\ ). 
The proof follows by writing 

/(V) -D =/(V) -D + (D - D) =/(V) - /(V) +(D - D) 

n n 

= ^P^^ - P(j)if(^^) - /(Aj))Vi ^ V,- D + ^ [A./'(A,) -1]Y,M V, D. (4-7) 



i,j=l i=l 



Hence, 



I /(V) -Dp = ^ [(p,, - 4) (/(A.) - /(A,)) tr(V,DV,)]' + ^ [(A,/'(A,) - l) tr(V,D)]^ 

ij'=l i=l 

where the two sums on the right hand side are the corresponding terms in (l4-6p . 
Therefore 

Lemma 4. For every Eulerian strain measure f there is a unique spin which minimizes the 

o 

difference between /(V) and D, and that spin is . The minimal difference is 



/(V) -Dp = J2 [(A./'(A.) - 1) tr(V,D)]^ (4- 

1=1 
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The proof follows using Lemma [3] in the form 

I /(V) -Dp > |D - D|2, (4-9) 

with equality iff = . 

4B. The logarithmic spin. Lemma [4] implies that the corotational rate of strain equals D if 
the strain measure has the property xf'{x) — 1 = 0. The only solution satisfying the condition 
/(I) = is f{x) = Inx, and the associated spin follows from (031) as = W + P'°sd where 

>g = I + (4-10) 

In A, -In A, Af - A| ' ^ ^ 

ri'°s is the well known logarithmic spin [19j. Hence, of all possible rates and of all possible 
Eulerian strain measures only the combination of the Hencky strain and the rate defined by the 
logarithmic spin together yield the strain rate D. This is the unique relationship between In V, 
D and which makes both the Hencky strain and the logarithmic spin special. This result 
was first derived by Xiao et al. [19] and may be summarized as: 



Lemma 5. The strain rate D is recovered only as the corotational rate of the Eulerian strain 
eC^) = In V with spin where the fourth order projection tensor P'°^ is given by eq. 1^-5] ) 
with / = In. That is, 

°log 

lnV=D. (4-11) 

4B.1. Some properties of the logarithmic spin. An instructive alternative form for P'°s is ob- 
tained by introducing 

pi° = (VKII-IK V)*Aln(V), (4-12) 

so that 

plog ^ pin + p£; = pin + p^ + p-R. (4.13) 

Each of the projection tensors may be expressed in terms of matrix elements pij = —pji according 
to eq. ([3^ : 

Aj - A,- n; Af+A^ , 2AiAj _ 1 



X^ + Xj' X^-Xy A2-Af In A, -In A/ 

The form of agrees with the formula for P^°^ derived by Xiao et al. [19] eq. (41)]. Note 
that 



sgnpj/ = sgnpjj- = - sgn(p.^.*^ - p^^). 

The implications are twofold. The first equalities indicate that the spin induced by both fi'°^ 
and by fi^ are in the same sense, relative to the underlying spin W. The second equalities 
imply that the relative spin induced by 0^°^ is of smaller magnitude than that of fi^. 

The f-spin, which is uniquely defined by the strain measure /, defines the skew matrix elements 
p\-. Consider the reverse problem: given some objective corotational rate defined by elements 

Pij, is there a function / such that p\- = pijl There is no such function for the spins fi^, il^ 

and ri'^, as the reader can readily verify. Obviously, / = In for Q, = 0^°^, but it remains an 
open question for general whether a strain measure function exists such that = Vl^ . 
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5. EULERIAN CONJUGATE STRESS-STRAIN PAIRS 

5A. Arbitrary strain and corotational rate. It was noted in the Introduction tliat tlie 
concept of work-conjugate stress-strain pairs is more complicated for Eulerian quantities owing 
to the fact that the connection between the strain rate and the stretching tensor is not evident 
a priori. This issue was resolved by Lehmann and Liang [8] who introduced the notion that the 
Eulerian pair r and e are defined to be conjugate if 



zi; = tr (QrQ* QeQ*) , (5-1) 

for some rotation Q. This clearly generalizes the Lagrangian work-conjugacy condition p-2|) . 
but it is necessary because of the fact that Eulerian rates are not as restricted. The definition 
(|5-1|) is equivalent to 

w = tr(T e), (5-2) 

where e= e + {e,fi} and = — Q*Q. Equation (j5-2p is taken as the starting point since it 
depends only on the corotational rate through the spin fi, i.e. Q is not required. 

For a given strain measure e = /(V) and corotational rate = W + PD the strain rate e 
follows from Lemma 13-31 The stress r is therefore conjugate to e if the following holds for all 
stretching tensors D, 

tr(r[V/(V)] Q(V)I)) =tr (crD). 

The fourth order tensor V/(V) is invertible for all strain measures. The necessary and sufficient 
condition required to determine r is therefore that the fourth order tensor Q is invertible. This 
requirement was obtained by Xiao et al. [T9] in a slightly different manner; basically that the 
six elements qij of eq. (j3-10p are all non-zero. Hence, are bounded, and exists. We 
refer the reader to Xiao et al. [19] for further details. 
In summary, the conjugate stress is 

T = [A/(V)] Q-V, 

where the order of [A/(V)] and is arbitrary since they commute. This is Theorem [H 

5B. Conjugate stress and the f-rate. An alternative approach is suggested by eq. (14- Th . Let 

e= FD, then the fourth order tensor F is by assumption invertible and the conjugate stress is 
simply T = F~^<T. The tensor F can be obtained directly in spectral form from eq. (j4-7p and 
easily inverted, to give: 

Lemma 6. For arbitrary strain measure and rate the conjugate stress can he expressed 



n 

y Vi M V,- cr. (5-3) 

.5^il+[fe.-4)(/(A.)-/(A,))]-' 
«7^i 



The conjugate stress satisfies 

|r - crp = |t - cr^P + |cr^ - crp, (5-4) 
where the modified stress tensor is 



n 

The proof follows from eq. (j5-3p by analogy with the proof of Lemma [31 Hence, |t — crp > 
\af — <Tp with equality iff = fi-^, and we deduce: 
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Lemma 7. For every Eulerian strain measure /(V) the corotational rate of the f-spin 
minimizes the difference between the conjugate stress and the Cauchy stress. The conjugate 
stress is then t = and the minimal difference is 

n 

l-H'-|:[(^-i)Mv..)]^. 

This proves Theorem ([2]). 

In general cr^ is not equal to the Cauchy stress for any strain measure, with the exception 
of / = In, discussed below. It is however, possible for af and a to coincide under special 
circumstances: if the three elements tr(Vjcr) simultaneously vanish. This is by definition a state 
of pure shear [11]. Hence, we have 

Lemma 8. // the Cauchy stress is a state of pure shear with diagcr = in the principal axes of 
V, then 

= a. 

The stress conjugate to /(V) equals the Cauchy stress if the f-rate is used. 

If the material is isotropic then the stress and strain share the same triad of principal axes. 
In that case diag cr expressed in the principal axes of V is simply the principle stresses, which 
vanishes only in the absence of stress. Hence the circumstances under which Lemma [8] applies 
cannot occur for isotropic materials. If the material is not isotropic but we restrict attention to 
linear anisotropic elasticity, then diag cr expressed in the principal axes of strain e will vanish 
only if both stress and strain are zero. This follows from the assumed positive definite property 
of the strain energy, equal to ^tr((Te). In summary, the circumstances under which Lemma [8] 
apply require nonlinear and anisotropic elasticity. This does not eliminate its possibility but it 
makes it difficult to envisage a situation when Lemma [8] would occur. 

5C. Logarithmic rate. The logarithmic rate, as noted before, is a special case of the f-rate. We 
conclude by examining the conjugate stress for arbitrary strain measure using the logarithmic 
rate. Xiao et al. ^19j showed that the logarithmic rate is the only one with the property of 
Lemma O i.e. among all strains and all rates, only In V and correspond to the stretching 
tensor D. This fundamental result for InV is generalized to arbitrary Eulerian strain measure 
e = /(V) by 

(V/(V)) (AlnV)D, 

which follows from Lemmas [2] and [H Now require that the work-conjugacy identity tr(T e^) = 
tr(<TD) holds for all D, and use the invertibility of the fourth order tensors V/(V) and AlnV 
plus the property that they commute. This implies that the stress conjugate to the Eulerian 
strain e = /(V) is 

T = ( A/(V)) (V In V) a for = 
It is straightforward to show that this can be expressed in spectral form as 

r = + X: ( : - 1) V. ^ V, a for = fi^"^. (5-5) 

This identity, although valid only for the logarithmic rate, shows how the conjugate stress in 
that case is related to the modified stress cr-f . The latter depends upon the strain measure 
/, and is optimal in the sense of best for all possible strain rates. Equation ()5-5p shows that 
the logarithmic rate is not optimal since r satisfies eq. (j5-4p with both terms on the RHS of 
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the latter non-zero. However, when the strain measure / reduces to hi then ^ a and the 
sum in (15-5P vanishes. This again shows the combined properties of the Hencky strain and the 
logarithmic rate as being doubly optimal for all strain measures and spins. 

6. Conclusion 

We have examined the implications of work-conjugacy with emphasis on Eulerian stress-strain 
pairs. There is, however, remarkable similarity in the form of the dual conjugate stresses for 
both Lagrangian and Eulerian strains. The similarity is evident from the identical format of eqs. 
(jl-5p and (|l-10p . which involve fundamental stresses T*-^^ and r^^^ defined by the strains /(U) 
and /(V), respectively. The Lagrangian stress T^"*^^ is called Biot stress or Jaumann stress, but 
there does not appear to be a common term for its Eulerian counterpart t^^\ 

The major distinction between Lagrangian and Eulerian work-conjugacy is that the latter 
requires the introduction of the corotational rate, which itself is quite arbitrary. We have shown 
that every permissible Eulerian strain measure /(V) has associated with it a unique corotational 
rate, the f-rate. The conjugate stress obtained using the f-rate is optimal in the sense that it 
is the closest possible to the Cauchy stress <t. The optimal stress, , is defined by / and a 
through eq. (j5-5p . and it reduces to the Cauchy stress if and only if / = In. This reinforces the 
results of Xiao et al. [20] for the the logarithmic rate and the Hencky strain, while generalizing 
the notion of the logarithmic rate to arbitrary strain functions through the strain dependent 
spin ri-^. 
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Appendix A. The spins Q.^, Q.^ and 9.^ 
From the definition of Q,^ in ()3-6p . and using F = RU, we have 

L = FF~^ = + RUU^R*. 

The symmetric and skew parts of this relation yield [16] 

D = ^R(UU"^ + U"^U)R*, W = + ^R(UU"^ - U^^U)R*. (A-1) 

Equation ()A-ip i may be solved for U in the form given by eq. (jl-4p . Substituting U in (jA-ip ? 
gives 

n^ = w-F(iKv-vKii)(iKv-FVK iy^T>. 

Let Vj, i = 1, . . . , n < 3, be the principal axes of B and V. The twirl Q,^ defines the rate of 
rotation of this triad by Vj = fi^Vj. The rate of change of the eigentensors of B follows from 
Vj = Vj (g) Vj as Vj = {f2^, Vj}. Equation (IS-ip o then gives 

n 

B = J^/3jVj + {n^,B}, (A-2) 

i=l 



16 



ANDREW N. NORMS 



which can be considered as an equation for Q,^ , similar to eq. p-3p . The solution follows from 
eqs. (|2^ and (f3^ as 

n 

= (I K B - B K I)* (B - ^ AV,) (A-3) 

i=l 

= W+(IK1B-BKI)*(IKB + BK1 1)D. 

Hence fj^ = W + P^D where is given by ()3-7bp . The rate of change of the principal 
stretches are obtained by substituting Cl^ into (IA-2p . as 

n 

= [ I - (I K B - B K I)*(I K B - B K I)] (I K B + B H I)D. 

i=l 

Then using eqs. (l2-6p and ()2-8p . we obtain the well known result 

n n 

^AV, = 2^AViKV,D ^ = A, tr(V,D). 

i=l i=l 

The twirl tl^ defines the rate of rotation of the Lagrangian principal axes Uj, z = 1, . . . , n as 

Uj = ri Uj. Hence, Uj = {ft ,Ui} where Uj = Uj Uj are the eigentensors of U. Taking the 
rate of change of the identity Vj = (RKlR)Uj linking the Eulerian and Lagrangian eigentensors, 
gives 

U» = (R K R)-^ {ft^ - n^, Yi}. 
The Lagrangian twirl is therefore 

= (R K R)"^ (fi^ - Q,^) = (I K - K I)* 2(U K U) (R Kl R)"^ D. 

This is related to the spin = W+ P^D defined via P^ of eq. (I3^ by = (RKR)"i(ri^- 
W). 
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